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ABSTRACT 


Velocity and temperature fluctuations from multi-level 
measurements in the marine boundary layer were examined to 
‘determine the validity of present formulations and prediction 
techniques, Spectral distributions were analyzed using an 
analog analyzer as opposed to the digital computer methods 
previously utilized at the Naval Postgraduate School. 

In general, the height variations of the dissipation 
rate of turbulent kinetic energy, momentum flux, and the 
temperature structure parameter observed in this study agree 
with the overland expressions developed by Wyngaard. Most 
deviations from expected resultS were attriouted to apparent 
anomalies in spectral results. 

It has been shown that the analog spectral analysis 
performed in this study produced information which is as 


reliable as that from digital methods, 
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I, INTRODUCTION 


This is a study of velocity and temperature fluctuations 
in the marine boundary layer. The data examined were obtained 
from multi-level measurements with tungsten wire velocity 
probes and platinum wire temperature probes, The measure- 
ments were made from the Naval Postgraduate School research 
vessel, ACANIA, Mean wind and temperature measurements are 
also utilized to complement and analyze the fluctuation data, 
A separate investigation on the mean data has been completed 
by Smedley (1975). 

The purpose of this study is to examine properties of 
small scale temperature and velocity fiuctuations in the 
fietne CHhys ronment, On the=™basis of spectral distributions. 

A Significant aspect of the study is that spectral distribu- 
tions were obtained utilizing an analog spectrum analyzer, 
rather than digital computer methods. The latter has been 

the method used previously in atmospheric turbulence experi- 
ments at the Naval Postgraduate School. The study includes 
descriptions of the variation with height of (1) €« (the viscous 
dissipation of turbulent kinetic energy) and (2) Cr (the 


temperature structure parameter), 


Le: 





II, THEORETICAL CONSIDERATIONS 


A. BACKGROUND 

This study is one of a continuing investigation of 
multiple-level overwater observations from a ship. The ship 
is the R/V ACANIA, shown in Figure 1, operated by the Depart- 
ment of Oceanography, Naval Postgraduate School. Because 
these were shipboard measurements, only high frequency small 
scale fluctuations could be measured, The following paragraphs 
will review the theoretical considerations used to interpret 
the small fluctuation data, Parallel reviews of this theory 
have been provided by Johnston (1974) and Bone (1974). The 
latter involved initial examinations of data obtained with 


the same instrumentation used in this study. 


B. TURBULENT KINETIC ENERGY BALANCE EXPRESSION 
In a simplified form, the turbulent kinetic energy 


balance expression for neutral conditions can be written 





— = € 


where u represents the mean wind at height z and the term 


su is the wind shear, Equation (1) states that mechanical 


production of turbulent kinetic energy, -u'w' i equals the 
molecular dissipation of turbulent kenetic energy, €. 
Equation (1) is often written as 
vo 28 oe (2) 


OG 


afl 








The R/V ACANIA 


Figure 1. 
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where uz = -u'w', and -u'w' represents. the vertical transfer 
of horizontal momentum by turbulent fluctuations. 
For neutral conditions in a constant stress boundary layer, 


the following expression is valid for the wind gradient 
IZ KZ 


where K is Von Karman'ts constant. Combining Equations (2) 


and (3) yields an expression to estimate UW, on the basis of 


* 


measurements of € at a known height Z. 


a 
Ue Soke (4) 
Equation (4) can be rewritten 
u 
* 
Ln € = -Ln z + Ln x (5) 


which is a form used in later analysis of the variation of 
—€ witn heignt, 
For non-neutral conditions the turbulent kinetic energy 


balance expression (Garratt 1972) is 


2 di 
ils CGR a 


For an unstable case CR. 90), 


6 (6) = (1 - 16 &)74 (7) 
and §€ = R.. 
For the stable case (R, > 0), 
bE) = 1+ 5E (8) 
and 
R 
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where Ri is the Richardson number, | 


= a ori 2 
R, = e/@ (3T/,,) / (30 /,,) 
Equation (5) can be rewritten as 


3 
— ee * 2D) . 
eee Ga) -6))* 29 5 


A plot of Ln ¢€ versus the first term on the right hand 
Side of Equations (5) or (9) would yield a straight line with 
-l slope. Figure 2 illustrates the possible effects of 
instability on the predicted plot of Ln z versus Ln e€, where 
Ln >, (R;) is an empirical function of the Richardson number, 
Deviations from this relationship couid be interpreted as 
either due to deviations from the constant flux assumption 
arising from the influence of wind-wave coupling or non- 
neutral stratifications, 


From the above discussion, ¢ can be used to estimate uy 


values... ie is an important parameter in boundary layer studies 


Since it defines the intensity of turbulent transfer of 
momentum toward the surface. It is assumed to be constant 


with height in the first 30 meters. 


C. SPECTRAL SIMILARITY THEORY 

Kolmogorov's second hypothesis is the basis for estimating 
€¢ and Cr (the temperature structure function parameter) 
from spectra of velocity and temperature fluctuations. It 
is based on similarity expressions for the inertial subrange, 
which is the wave number region where energy is transferred 


without dissipation from lower to higher wave numbers. 
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Fieure 2, An Iliwstration of the Effect of Stability on 
Vattiacion of ©€ with Height, 
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The following relationships were proposed by Kolmogorov 


for spectral densities in the inertial subrange, 
_ Dy | 
OG —<, « kK , (10) 


on(k) = ¢, C2 ke (11) 


where $(K) and On (Kk) are the one dimenensional spectra for 
velocity and temperature respectively and k is the streamwise 
component of wave number as defined below. 

€ has already been discussed in it's relationship tou,. 
Another turbulence parameter of interest in this study is the 


temperature-structure parameter, C2, which appears in 


TT 
equation (11). It is important because it is related to the 
refractive-index-structure parameter, oe. by the expression 
Ae as ™~ 6 2 
c2 = “a my oO D C2 (12) 
12 T : . 


Cx is a primary parameter in the characterization of optical 
propagation in the environment, 

Height variations of Cr are dependent on stability, as 
was the case for e«. Predicted height variations for various 
stability conditions are shown in Figure 3. Wyngaard, et. al. 
(1971), noted that Cr decreased with height as g7*t* sor 
stable conditions and as Z~*/? for unstable conditions. 

Since velocity and temperature fluctuations are measured 
at a fixed point in the flow, the resultant spectra are 
realized at a "temporal" frequency, f, To obtain an ¢e and 


Cr the time (f) and.space (k) scales must be related by 


Laviomuem Gloss) irozen turbulence’ hypothesis, 1,e€., 
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Figure 3. An illustration of the effect of stability on 


the variation of Cr with height. 
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where u is the mean wind speed at the measurement level. The 
term "frozen" implies the turbulence pattern remains unchanged 
as it sweeps past the probe. 


Equations (11) and (12) can also be written as 
fo(f) = ko(k) = c,e?/*K?/3 | (13) 


fop(f) = Kbp(K) = 0, c2K72/3 (14) 
where o(f) and Op (t) are spectral density values with units 
of (m/sec)?/Hz and °C?/Hz respectively, 

Final forms of Equations (13) and (14), in terms of 9$(f) 


‘and On(f), which are measured, are 





/ 2/3 Nesi/e2 
s = \° a) (19) 
/ ae k2/3 
ee il f(t) (16) 
2 
PRere  cipimtcaliy, c,= 0.5 and c, = 0.25. 


From Equations (15) and (16), with measured f, $(f), d(T), 


2 


Jk 


These were the final forms used to determine the variation 


and u, € and C, can be determined at each level of interest. 


with height of ¢€ and oF values. 
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III, DATA COLLECTION AND SENSOR CALIBRATION 


A, THE PLATFORM 

Observations were obtained aboard the R/V ACANIA anchored 
off Monterey, California in Monterey Bay, Measurements were 
made at four levels on two masts spatially mounted on the 
forward deck of the ship. Figure 4 illustrates the sensor 
locations on these masts. Figure 5 shows the vane and probe 
arrangement, The vane maintained a one jamemetiomel profile 


of the wind during measurements, 


B. INSTRUMENTATION 

The velocity fluctuation sensors were TSI model 1210 
probes with sliding support shields and tungsten wires, The 
shields permitted isolation of the sensing region for deter- 
mination of the undisturbed velocity, Se before and after 
each experiment. The platinum coated tungsten wire was small 
enough to resolve the viscous dissipation scale without any 
need to apply wire length corrections. The wire was fitted 
with plated ends for isolating the sensing area and thereby 
minimizing flow disturbance. 

Electronics associated with each probe were a TSI Model 
1054B Linearized Anemometer, a TSI Model 1056 Variable Decade 
Module and a Model 1057 Signal Conditioner. The signal 
conditioner was only used during the March 1974 observations. 


The anemometer had a linear frequency response from DC to 10 


1138) 
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Sea Surface. 


Figure 4. Mounting arrangements, 
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Vane and probe arrangement. 


Figure 09. 





KHZ and the variabie decade module operated with a O-60 ohm 
range. The signal conditioner had a linear frequency response 
from DC to 400 KHz and permitted voltage suppression in one 
volt steps from 0-29 volts. 

sensor placement required exceptionally long cables, but 
any decrease in system frequency response had little effect 
in the frequency band of interest, 

In general the performance of the hot-wire system was 
excellent throughout the experiments, 

Temperature fluctuations were also measured with a TSI 
Model 1210 configuration using platinum wire sensors. 
Electronics used were the GTE Sylvania Model 140 thermosonde 
temperature modules, which can be used to make both single 
point temperature fiuctuation and AT measurements. The 
temperature module is designed for measurements of temperature 
variations down to 0,001°C at frequencies up to Ee 

Mean measurement systems are described in smedley's 


GLO7S) study. 


C. SENSOR CALIBRATION 
1. Hot Wires Calibration 
In-situ calibrations of the velocity sensors were 
accomplished several times during each experiment, In this 
procedure, recordings were made of both the cup anemometer 
RPM count and the corresponding hot wire voltage output. The 


sensor wind speed and voltage relationship is given by 


_ 2 
> (CT See C7) 


22 


where v is the PenLeaGe and u is mean wind speed for a partic- 
ular level, The constants a and b are calibration curve slope 
and intercept respectively obtained from the in-situ calibra- 
tion. Figure 6 is a sample calibration curve constructed 
from in-situ values in one period, 

To convert the power spectral density levels to 


velocity Units required a calibration factor such that 


Cay =u" (18) 


: . : ; em/sec 
where cis the calibration factor in cm/sec 


a 
volt» Vv is the 


voltage flucuation and u' is the velocity fluctuation, 


Differentiating Equation (17), we get 


ay 
(44 ) or (19) 


or 


Cape (20) 


Toeoneain the calibration factor, c, Equation (17) 
was initially solved for v using the appropriate sensor 
level's mean wind speed u and also the particular sensor's 
calibration curve slope and intercept. 

2. Platinum Wires Calibration 

Wire resistance is related to temperature by the 
equation 

Reel adaeo( f= TO) ) (2a: 
where R = resistance, T = calibration temperature, To = O2e- 
Ro = resistance at To, and a = .0039/°C (temperature 


Coctitvervent Of resistance for platinum). Each wire’s 
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Figure 6. 


Hot wire calibration plot, 
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resistance was measured in the laboratory at a temperature T 


? 


to determine Ro, 


We can relate temperature fluctuations to wire 


resistance fluctuations by differentiating Equation (21), 


dR Geto al 


or 


Rio- a ko T'! (22) 


Recorded voltage fluctuations can then be related to 
resistance changes by the transfer function, H, in the 
expression 

VR” (23) 
where H was determined for the four bridges by varying the 
resistance of a TSI 1056 Variable Decade Probe Resistance 
and recording the corresponding voltage changes, Figure 7 
is a calibration curve of voltage versus resistance, where 
the slope determines the transfer function, H. 

Finally, to relate recorded voltage fluctuations to 
temperature fluctuations substitute Equation (23) into 
neiac ton (22), Siving 

= (ea) . a 
Or 


me vi (25) 


where Ch is the temperature calibration factor in [Comte 
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Figure 7. Platinum wire transfer function plot, 
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D. SIGNAL PROCESSING 

The data BEOmet Mew Loum Sevelsowas recorded simuleancously, 
in analog form, on a 14 channel tape recorder, 

The signals were continuously monitored with an 
oscilloscope and an 8 channel strip chart recorder to insure 
the quality of the recorded signal and to detect equipment 


failures. 
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IV, DATA PROCESSING AND ANALYSIS 


A, PRELIMINARY ANALYSIS 

Selection of data periods was based on signal evaluations 
of strip chart records, mentioned in section III.D, and on 
evaluations of experimental logs of equipment status and 
ship positioning. Primary consideration was given initially 
to obtaining 18 minute records with temperature and velocity 
Signals at all levels, This was later chaneed cem2l minute 


records due to characteristics of the anolog spectrum analyzer, 


B. PROCEDURES 

Several procedures were adopted to best handle the 
turbulence data. These are outlined in Appendix A in detail. 
All spectral analysis was performed by processing the recorded 
analog signal on a Federal Scientific Model UA-500-1 Spectrum 
Analyzer-Averager shown in Figure 8. 

The averager output was a logarithmic plot of the power 
spectral density in Volts*/Hz, versus frequency f in Hz. The 
plotter was calibrated so that all plots would concide with 
Bone's (1974) digital plots. Thus, a preliminary visual 
comparison was made to verify the spectra prior to computations, 
Figure 9 illustrates a sample of the plot comparison for a 
data period common to Bone's (1974) study and this study. 
Figure 10 shows a typically good graphical plot of spectrum 
versus frequency used in this study. Figure 11 is an example 


of a poor spectral plot. 
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LOG, FREQUENCY 


Comparison of digital and spectrum analyzer 
velocity spectra for 1855-1913L, 18 Jan. 74. 
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C, SPECTRAL INTERPRETATION 

scaling of the spectral plot to relate an RMS input 
voltage to a power per unit of frequency is discussed in 
Appendix B, | 

The single point method as previously discussed by 
Johnston (1974) and Bone (1974), was used to obtain ¢€ and 
Cr values from the velocity and temperature spectra. Using 
many spectral points in the inertial subrange (-5/3 slope 
region), an intercept at an arbitrary frequency f can be 
obtained. This intercept value yields the power spectral 
density (PSD), after correction. 


1. Determination of ¢« and Wels 


The measured PSD value was converted to a spectral 


aGensity in engineering units by the equation 


c? + PSD 


o(f) 


(cm/sec/ . OES = (em/sec)*/,,, (26) 


2 
volt? 
where c is the hot-wire calibration factor, The solution to 
Equation (15) yielded a value for ¢ at each level. Then, 
Equation (4) and (6) were solved for a u, value and a 


corrected U, value respectively. Corrected U, values were 


*K 
not valid for Richardson numbers > 0,2. 


2. Determination of Cr 


The following equation was used to convert PSD 


values from measured spectra to dimensional values 


on (f) CH + PSD 


O 2 2 —_ VWA2 
( Oi ross) = Voit ly ESS ie © 
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where Ch is the calibration factor for the temperature 


system. Solving Equation (16) for Cr yielded a value for 


each level Z, These values were then applied to Equation 


(12) to obtain Cy for possible comparison with optically 


derived Cy values, 


34 





V, RESULTS 


A. DATA PERIODS EXAMINED 

Mmnulcthliatane velocity, uw’, and temperature, T’, data were 
studied for the nine periods shown in Table I, Four levels 
of both velocity and temperature were desired, but temperature 
Signals were not always available for all levels due to 
various component failures; Table I summarizes the periods 


examined and what data was available for each period, 


B. VELOCITY SPECTRAL RESULTS 

For seven of nine periods studied, four levels of 
velocity, u', were available. Equation (15) was used to 
calculate ¢ from values at each level from selected frequencies, 
f, spectral densities, $o(f), and wind speed, u, Each computed 
value € was then used in Equations (4) and (6) to determine 
vel Homerae  irict1on velocity, u,, 

Table II summarizes the e« and u, results along with 
stability corrected values of U,, profile u, estimates 
obtained by Smedley (1975), and Richardson numbers, Ra. Most 
spectra exhibited the -5/3 slopes for the inertial subrange 
region which validates the method used to determine e from 
the spectral densities, Except for one value on March 27, 
Weynvatwes are alll less than 0.0 Wess ” and most are in the 
expected range of 0.1 - 0.5. Most periods of u, seem to 


reflect the constant flux theory accurately, especially if 
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TABLE II 


€ and ge Results 


HEIGHT C vt Ae 7 R. 


= 5 = * *K 
(m) (m-—sec °) (m-sec /) (cCorrecrved) (profile) z 

Period 1 Zoe warch 74 1841-1902L 
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FOS .019 ,38 28 .16 
13.8 Oa 4b * Rests) 
Period 3 Zfeaktarch 74 1928-1945L 

Aan . 240 2 ,65 .26 , O06 
6.6 .039 .46 eal abs 
7.6 ,018 aot eis. wo 
ILS) as 210) 2 31 48 * yoo 
Period 4 27 March 74 2117-2138L 

4,1 044 ool 34 cio, a 09 
6.6 025 ,40 0 . O06 
35 .016 ANA! LONG. vey 
Period 9 28 March 74 1902-1923L 

AL al 039 76) .38 610) Oz 
Cae 0030 .20 a le’ 09 
7-6 0030 5A 0) .19 .05 
13.8 .0013 KS: plece pails 
NOTE 


: computed from equation 19 


computed from equation 6 


* no data 
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TABEE Tie comtpanued) 


u R 


uy, u : 


HEIGHT Ee 
Period 6 

4.1 ,020 
B26 O13 
7.6 ,010 
3.8 O11 
Period 7 

4 dl .0089 
6.6 ~0015 
736 ~00049 
Ley ete! 019 
Period 8 

abl .0099 
6.6 - O00G3 
7.6 ~0017 
iS. 5 ~0039 
Period 9 

a) 0077 
6.6 6076 
720 0013 
io. 6 ,045 


de 
CCOrte) 
2059-2117L 
PRS eau 
woo 14 
Gol : * 
.38 * 


18 September 74 1643-1701L 


~29 26 
29 26 
re erllera 
27 .26 


23 24 
27 28 
ne, 16 
»62 64 
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Obviously bad levels are discarded in the.periods, Those 


results which prevent nearly the same u, flux for all 


* 
levels in a period, were readily attributed to a departure 
neon tne =-o/3 spectral slope. A comparison of this study's 

u, and stability corrected u, values with Smedley's regression 


u, values, provides close agreement for at least two levels 
for each period, The comparatively high u, Values fomplevel 
4 in the stable cases has not been explained, but it also 
occurred in one of the two unstable cases. 

Pieunre la 92S a typical plot of height versus ¢« for a 
Stable case, It shows close agreement with the -1 slope 
expected for the near neutral case, according to Equation (5), 
Figure 13 illustrates the unstable case. Here again we have 
agreement with tne neutral case -1 slope, except for level 
3. Poor results for level 3, in this case, can be related 
to a poor spectral plot with insufficient amplitude and a 
slight departure from a -5/3 slope. This deviation is not 
expected since level 2 and 3 are within one meter in height. 

Composites of height versus ¢€ for all periods are shown 
in Figure 14 for stable and unstable cases. Obviously bad 
data was not included in these composites. Both the stable 
and unstable cases show close agreement to the near neutral 
prediction, This differs from the previous results of 


Bone (1974), in that his slopes deviated from the -1l value, 


as expected for non-neutral conditions. 
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C. TEMPERATURE SPECTRAL RESULTS 

Ciemuetrod had four levels of temperature, T', data, 
while eight periods only had three levels of T' data. Cr 
values were computed by applying Equation (16) to values at 
each level for selected frequencies, f, spectral densities, 
on(£), and wind speeds, u. Equation (12) was then applied 
to the resultant Cr tomeet Co. The Cx values will be used 
by investigations into coincident optical propogation 
experiments, 
i and C2 results for this 
study. It also includes Smedley's (1975) profile values of 


Table III is a summary of the C 


Cr for a comparison, The height variations of Ca rye eae 
stable periods were generally consistent with predictions. 
During unstabie periods there appeared to be some iacK of 
consistency between levels. Here again, most inconsistencies 
were traced back to a poor spectra. Cr values which are 
within an order of magnitude of Smedley's Peels are consid- 
Bred a fair=correlation. Such correlation does not occur for 
Gipeebevol in peraeds 1, 4, and 7. 

A typical plot of height versus Cr for the stable case 
is shown in Figure 15. It agrees very well with the expected 
-3/2 slope, The unstable case is shown in Figure 16, The 
expected slope of -3/4 for height versus Cr agrees reasonably 
well with two of the three data points. The 4 meter level 


point can be attributed to a spectra with a slope greater 


than -5/3 along with comparatively small amplitudes. 
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TABLE III 


2 2 
Cr and Cn Results 
HEIGHT oom Gre on R, 
(m) (Sc? m7?/?) (profile) (m7? /?x1071*) 

Period 1 2% March 74 1841-1902L 

4,1 .0017 * ,16 ~25 

6.6 0034 > woe OG 

7.6 ~0034 *K nom ,62 
Period 2 27 March 74 1902-1923L 

4.1 ~0021 7,045 20 .06 
6.6 ~0022 ~O061 es i4 

7,6 .0043 .0032 nae 16 
Period 2a Mareh 74 1928-1945L 

Al ,00a7 .042 55 .06 
626 ~0050 ~0065 .48 et 

7.6 20062 ,0036 .60 mes: 
13.8 ZOOLS * wlll SS: 
Period 27 March 74 211s ook 

ak .0036 OFZ 0 O09 

6.6 BOOS 2 .020 oa l .06 
8) are. Fools .00089 17 17 
Period Zor March 74 nOO2— lo 2 ois 

4.1  OIaL % ,016 ta Fe 202 

ES POzG .0040 2,0 Poo 
13.8 .0018 ,00055 .18 ~13 


A 


7 value computed from Equation (16) 


- values from Smedley (1975) 


K 
profile unavailable, R. >, 20 
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HEIGHT 


Period 
4.19 
7,64 

13.89 


Period 


TABLE III (continued ) 


2 2 
Cr Cn 
(profile) 
March 74 2059-2117L 
.014 SA, 
. 0046 3.0 
* 2nee 
september 74 1643-1701L 
»0016 029 
~0012 015 
5 013 


september 74 


.0063 
FOO 
.O015 


September 74 
,O08l 


. 0036 
.0016 


1735-1753L 


. 048. 


24 
,O7 


1900-1918L 


D 
we 


ea tS Te 
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Figure £7 includes composites of height versus Cr 1eOTe 
al serious for both the stable and unstable cases. For 
these composites, obviously unrepresentative data were not 
plotted. In the stable case, agreement with an expected 
Slope greater than -3/2 is possible, The lowest level shows 
the greatest range in relation to the slope line. This is 
attributed to a comparatively weak amplitude spectra in 
many cases. Actually, this level should have the largest 
amplitude spectra in comparison to the other three levels, 


There is good agreement with the expected -3/4 slope for the 


unstable case. 
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VI, CONCLUSIONS AND RECOMMENDATIONS 


A, CONCLUSIONS 

In general, the height variations of ce, U, and Cr 
observed in this study agree with the overland expressions 
developed by Wyngaard, The agreement with the constant flux 
assumption iS an important result of this study. Most devia- 
tions from expected results were able to be attributed to 
apparent anomalies in spectral results. Correlation of Ce 
results with Richardson number derived Cr values was not good, 
However, Since spectral results were generally consistent 
between levels, the spectral results are considered to be 
more reliabie. 

It has been athena that the analog spectral analysis 
performed in this study produced information which 1S as 
Feliable as that from digital methods. Application of the 
analog method on the large amount of available data should 


furnish information as to the reasons for identified 


departures from expected results, 


B, RECOMMENDATIONS 

Future experiments would be improved by relocating level 
3 further up the mast, say to 10 meters, to obtain a greater 
separation between levels 2 and 3. 

Improved strip chart monitoring of u' and Tv Signals 


during experiments is needed. Improved monitoring of these 
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records during the experiment could reveal anomalous signal 
levels and corrections during the experiment could be 
accomplished. In general, a decrease with height of both 
u*’ and T' signals is expected. 

More frequent in-situ calibrations of the hot wire 
anemometers should be made to increase the accuracy of results, 
especially for experiments of long duration. Frequent resis- 
tance measurements, on the platinum wires during the experi- 
ments is also necessary for this purpose, 

Finally, it is recommended that the spectrum analyzer's 
"vernier" input amplitude control be tested on a few poor 
periods to determine if the analyzer sensitivity could be 
part of the problem in the lack of consistency between some 
levels. Readjustment for each level and also obtaining a 
spectral amplitude calibration plot, would result in maximum 
sensitivity and minimum distortion of the spectral output, 
A.C. harmonics from the ships power supply should also be 


filtered out to prevent.spectral distortion, 
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APPENDIX A 


A. SPECTRUM ANALYZER-AVERAGER 
1. General 

The following discussion is for further amplification 
of the spectral analysis procedures used in reducing all data 
for this study and will serve to aid those who will continue 
to study turbulence data in the same manner. 

2. Equipment 

All spectrum analysis was performed by using a 
Federal Scientific Model UA-500-1 Spectrum Analyzer-~Averager, 

- which is a time-compression spectrum analyzer, using a digital 
speed up loon. 

In time-compression spectrum analyzers, the input 
data is digitized and stored in amemory. While the memory 
is being updated with new data, its contents are read out 
repeatedly at a very high rate, The speeded-up high frequency 
replica of the analog input signal is spectrum-analyzed, using 
a wide filter which is being stepped rapidly in real time. 

A recorder analog input signal is converted to digital 
after passing through an anti-aliasing low-pass filter, The 
sampling frequency is controlled by the analyzer frequency 
operating range. The digital signal is stored in a 1500 word 
digital memory, circulates in the memory, and can be non- 
destructively read out at an accelerated rate. This memory 


Output is then put through a digital to analog converter. 





The converter output is a time-compressed, accelerated form 
of the original analog input signal, It passes through a 
heterodyner and 1S applied’ to a stabllemerystal filter, causing 
the filter to observe successive frequency-resolution cells, 
A weighting function modifies the analyzer's frequency 
selectivity. Settling time of the crystal filter is short 
enough to obtain one spectral point each time the memory does 
one recirculation. Five hundred spectral points can be obtained 
in 50 milliseconds. Averaging of the spectral wave form is 
then accomplished with a digital ensemble averager to obtain 
statistical stability in mean spectral characteristics. The 
final result in storage is a real-time power-spectral- 
density (PSD) representation available at output terminals. 
The averager has two 500 word memories. 

A Hewlett-Packard X-Y Plotter was used to obtain a hard 
copy of results for each sampling period. 

3. Operating Procedures 

a. Analysis Range 
The range considered best for this study was 200 

Hz for all data analyzed. This resulted in a sample interval 
of 2.5 seconds (memory time). Since the highest frequency 
of interest was Peri: 100 Hz, the Nyquist sampling theorem 
was satisfied. The sampling rate for the above analysis 
range is 600 Hz, with a filter bandwidth of 0.4 Hz. 100 
percent of the input data is processed for the 200 Hz 


analysis range. 


93 





b,. Input Mode 
| Input signal coupling was used in the AC mode to 
analyze only the AC component of the input Signal, Hence, 
the analysis was not subject to DC offsets. This mode has 
a low-frequency 3-dB cutoff at less than 3 Hz, A presampling 
filter provides at least 50 dB of attenuation per octave 
fall off (i.e. 50 dB at 400 Hz). This filter and a sampling 
rate of three times the analysis range eliminate aliasing, 
c. Input Amplitude 
The input amplitude was adjusted for varying 
values of the RMS input voltage to attenuate the level of 
the input signal. The Vernier Cal. control was left in the 
detent position so that the input amplitude switch was 
functional from O-dB through -40-dB, in 10-dGB steps only, 
An RMS voltmeter was used to monitor sample portions of the 
input analog signal to determine a volts RMS full scale. 
The input amplitude was then adjusted for a volt setting so 
that full scale indications didn't exceed this volt setting, 
Periodic noise spikes were not considered. Overflow indica- 
tions were not normally excessive if this procedure was 
followed. This indicator may flash a few times per second 
without danger of significant signal distortion. 
dad, Transient Capture Controls 
Depression of the manual button, with the hold 
light off, enabled the memory to accept new data continuously. 
The two other available modes were not used in this analysis. 


theme the trigeer polarity and Sensitivity "conerols sana tne 
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transient auto delay switches are not functional. The 
external sampling switch must also be off. 
e. Output Spectrum Gain 
Pushbutton operation introduces a gain (O GB, 
+10 dB, + 20 dB) to the analyzer output prior to averaging of 
eecctra, Operation, for this study, was in the +10 dB (gain 
x 3.2) mode, due to the characteristics of the data signals. 
The output gain overload light indicates when spectral peaks 
are above full scale and are clipped, This occurred 
infrequently. 
f. Amplitude and Frequency LIN/LOG Switches 
Logarithmic scale for both amplitude and 
frequency were used for all spectra, 
g. Display Sweep Controls 
The FAST OSCILL mode must be used for the averaging 
of a signal and the SLOW PLOT mode should be selected for 
mlbecting the memory contents of the A=-Y pillotter. The START, 
STOP, RESET pushbuttons .control the averaging or plot mode 
elected. In slow plot an automatic rate control is used for 
large signal excursions to avoid plotter slew-rate problems. 
h. Averager 
(1) Active Memory Control, This switch selects 
that memory where averaging takes place and whose contents 
are displayed or plotted. The ERASE pushbutton will erase 
the selected memory if FAST OSCILL is selected, 
(2) <Averagine Method Controts..1n tne SUM miecac, 


each non-redundant analyzer spectrum 1S added to previous 


OOD 





spectra stored in the active memory, The INSTANT SPECT 
control can't be depressed when an averaging process is to 
be started. 

(3) Average Stop Mode Controls, The manual 
mode was used in this study to stop each analysis at 18 
minutes with the STOP control. The Number of Spectra-N 
oweecn then functioned as a gain control 1f set at N less 
than 256. This proved to be a problem, in that the analyzer 
looks at this switch to indicate how many samples are coming. 
If the number N selected and the number N or samples do not 
Beamcide when analysis stops, then the spectra is not 
complete. There is no problem in the Stop-At-N mode. 

ical bration 

Ali other operating procedures and the calivbration 

procedures of Analyzer-Averager and displays can be obtained 


from the manufacturer's instruction manual, 
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APPENDIX B 


A. SPECTRAL X-Y PLOT SCALING PROCEDURES 

To obtain a power spectral density level, corresponding 
to an RMS voltage input, a calibrated scale overlay had to 
be constructed for use in spectral interpretation. 

Amplitude scaling was accomplished by an internally 
generated Signal available in the AMPL, CAL input mode, 
When O0-dB (.iv) input gain was selected, a cal signal with 
amplitude approximately equal to 0,1 voits RMS appears on 
the display as a spike at 100 Hz indicating full scale 
amplitude. This spike was then plotted in the X-Y plotter, 
SueOne wich bench MameksS, uSinge a sSspeetral tain of Od (xi) 
Sabino oucCcesSively, the input gaan was sreppea to 
attenuate the amplitude of the spike by 10 dB, and a plot of 
the height of each spike was obtained on the X-Y plot. This 
plot now represented a graduation of RMS input from a minimum 
of 0.001 volts to 0.1 volts maximum. These voltages were 
then converted to LOG 10 and a graduated scale was constructed 
to that the LOG 10 of volts RMS could be interpolated from a 
Spectral et. The amplitude scale was adjusted for each 
Spectrum as a function of input attenuation and spectrum gain. 
These LOG 10 values were then converted to PSD levels for use 
tapecalenlating epsilon and Structure parameter values. me 


moderal Scientific Manual's equation to convery to Pavers 
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2 
0.79 (cal, level Veus) 


BoD elevonm@ly/ aa 
Filter Bandwidth 


where the filter bandwidth (8 = 0.4) is a function of the 


analysis range selected, 
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